Entanglement capability of self-inverse Hamiltonian evolution 
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We determine the entanglement capability of self-inverse Hamiltonian evolution, which reduces to 
the known result for Ising Hamiltonian, and identify optimal input states for yielding the maximal 
entanglement rate. We introduce the concept of the operator entanglement rate, and find that the 
maximal operator entanglement rate gives a lower bound on the entanglement capability of a general 
Hamiltonian. 
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Two-qubit unitary gates occupy a central role in quan- 
tum information science p| , particularly because of the 
capability of these gates to generate and enhance entan- 
glement of the state of the system. In general, entangle- 
ment capability can be enhanced by introducing ancillary 
states 0, El , but for the important and interesting case 
of the Ising Hamiltonian Rising = ®z ® <J Z it is ancilla- 
independent |j,|5j (where cr z =diag(l, — 1) is a Pauli ma- 
trix). The independence of the entanglement capability 
on ancillas is a consequence of the self-inverse property 

Rising = #ising- 

We generalize this result to all Hamiltonian evolution 
of the type 



H — X, 



X f 



(1) 



such that X % = X7 X g H t for i e {A, B} and H = H' 1 , 
with the Ising Hamiltonian being a special case. Here, 
we assume that Xi is not an identity operator 1,. Due 
the self-inverse property of the Hamiltonian, we have the 
evolution operator [h = 1) 

U(t) = e~ lHt = cost 1a® Is -isiat X A ®X B . (2) 

We employ operator entanglement 0,0 and operator en- 
tanglement rate to characterize the entanglement capa- 
bility of the self-inverse Hamiltonian evolution: these ap- 
proaches yield elegant analyses that are simple to apply 
for determining entanglement capabilities. We find that 
the maximal operator entanglement rate always gives a 
lower bound on the entanglement capability of a general 
Hamiltonian. 

The entanglement capability of a Hamiltonian 0, 0, 
is defined relative to a specified entanglement measure. 
We use von Neumann entropy as our entanglement mea- 
sure of a pure state | vf) S Hab ■ 



E{\y)) = -Tr A (p A \og PA ) 



(3) 



where pa=Ttb(pab), Pab = |*)(*|, and log is always 
base 2. The entanglement capability of Hamiltonian H 
is defined as the maximum entanglement rate when a 
pure state is acted on by the associated evolution op- 
erator U(t) = exp(-iHt). Mathematically, the ancilla- 
unassisted and ancilla-assistcd entanglement capabilities 



are defined as 0, Q 



Eh 

pane 



max|<p) EWAB r(t)|t^ > ( 4 ) 
su P|*)ew A , ABB , r Wlt-*o, (5) 



respectively, where T(t) = dE[U(t)\^(0)}]/dt is the state 
entanglement rate, and ancilla systems A' and B' are not 
acted on by Hamiltonian H . For the ancilla-assisted case 
the entanglement refers to the bipartite entanglement be- 
tween systems A' A and BB' . 

Let us first consider the unassisted case, and the ini- 
tial pure state of the system can always be Schmidt- 
decomposed as 



(6) 



l*(0)> = E 



where {|VVi)} and {!</>«)} are orthonormal sets of states, 
and A„ > Vn. As |^(0)) is normalized, J2 n ^ n = 1- 
The state at time t is described by the density opera- 
tor pab (i) = U(t)p AB (0)U^(t), which satisfies pab (t) = 
—i[H, PAB(t)]- Thus, the reduced density operator pA.(t) 
satisfies 

PA(t) = -iTr B [H, PA B(t)]. (7) 
From Eq. @ , we know that the entanglement rate is Q 



r(t) = -TY A [p A (t)io gPA (t)}. 



(8) 



Substituting Eq. (JJJ into Eq. ©, we obtain the entan- 
glement rate at t = as 

r(t)| t _ = iTr A {Tr B [H,p AB (0)]log PA (0)}, (9) 

which is general for arbitrary Hamiltonians and initial 
states. By varying initial states we can maximize the 
entanglement rate. 

Equation © is solved by first obtaining the results 



Ttb[H,Pab(0)] = E \J ^rn An (X B ) mn [Xa , \lpn)(^m\], 
ran 

(10) 

where (X B ) mn = {4>m\X B \4> n ) , and 

logpA(0) = ElogA n |^„)(V„|. (11) 
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The entanglement rate at t = for our Hamiltonian H 
is then given by, 

r(t)|t_>0 = V^m^n log T 21 PCt ) tnn (X B ) mn , (12) 

mn 

where (X A ) mn = (V> m |^4 |^n)- 

From Eq. (JT^J, we obtain an upper bound for the en- 
tanglement rate 



which 



= ^ ' (Am + 



log 



A TO 



|(X A ) m „| 



, A m 
logT— 

An 



X IPCOmnl \{X B )mn\ 

<7^(A m + A„) l(^A)m„| \(X B ) ■ru n | 
mn 

— 7 ^(^ m + (l(^)mn| 2 + |(Xs) mn | 



=/3, 



(13) 



where 



/3 =2 max y x(l — x) log[x/(l — x)] 
xe[o,i] 



=2V^o(l - z ) log[a;o/(l - aso)] « 1-9123, (14) 

and xq — 0.9128. The first line simply follows 
from the triangle inequality, the second and third lines 
from Ref. Q, and the fourth line from the inequality 
2|o6| < | a. | 2 + \b\ 2 for any complex numbers a and b. 
Finally, the last line results from the self-inverse prop- 
erties of Xa and Xb- To see this, let us examine the 

Sum J2mn Kn\{.XA)mn\ = J2mn ^m(XA)mn(XA)mn = 

Y, m > 'rnY, n {XA)mn{XA)nrn = 1, where the second 
equality results from the Hermitian property of Xa, and 
the last equality from the self-inverse property Xa — 
(Xa) -1 - We have seen that the self- inverse property is 
essential to obtain the upper bound (3. Another feature of 
Eq. ijT3|) is that the result is applicable for any pure state 
with or without ancillas. Therefore, the upper bound f3 is 
ancilla-independent. Next we show that the upper bound 
can be saturated by optimal input states. 

Since the self-inverse operators Xi, i £ {A,B} satisfy 
X? = 1, the eigenvalues of Xi are ±1 with the corre- 
sponding eigenstates (possibly degenerate) denoted by 
\±)j. Then, we construct the optimal input states given 

by m 



Mo)) 



(I- 



— xq 



\+)a-\~)a)®(\+)b-\-)b)- (15) 



This state is of Schmidt form with Schmidt number 2. 
From Eqs. JT2J| and IpHJ). we obtain T(t)\ t ^ = f3. There- 
fore, we can always find the optimal input states for 



E 



H 



0- 



(16) 



Result (|lf)H extends results for the Ising Hamiltonian 0, 
4] , and we do not restrict the dimension of Hilbert spaces 
on which the self-inverse operators Xa and Xb act. 
The Hilbert space can be finite-dimensional or infinite- 
dimensional. 

Self-inverse Hamiltonians not only exhibit elegant en- 
tanglement capability, but arc also physically meaningful. 
As an example of a physical system with optimal input 
state and maximal entanglement rate /3, we consider a 
spin-j system with a (2j + l)-dimensional space spanned 
by basis states {\n)j = \j;n — j), n = 0, 1, ... , 2j}. We 
define a number operator N = J z + j of the spin-j 



system satisfying Af\i 



Here, commutators 



for J z and ladder operators J± satisfy the su(2) algebra 
[J 2 , J±] = ±J±, [J_|_, J_] = 2J Z . Then, using the num- 
ber operators, we construct the self- inverse interaction 
Hamiltonian 



Hi = (-l)^O(-l) 



(17) 



where (— 1)^ is a parity operator of the system. This 
kind of Hamiltonian is realizable in physical systems • 
From Eq. i|17fl . we see that Hamiltonian Hi reduces to 
the Ising Hamiltonian -Hising for the case of j = 1/2. 
When the dimension d = 2j + 1 is even, i.e., j is a half 
integer, Hamiltonian Hi is equivalent to the Ising Hamil- 
tonian in the sense that the operator (—1)^ can be writ- 
ten as Ij/2 ® where Ijv denotes the N x N identity 
matrix. For integer j, Hamiltonian Hi is not equivalent 
to the Ising Hamiltonian. 

Let us consider the SU(2) entangled coherent 
state (ECS) 01 I 7 ?)ecs as an input state 



\v)ecs = V^olv) ® \v) + VI - sol - V) ® I - V), (18) 

where \rj) — exp(r?J+ — ri* J_)|0) 7 - is the SU(2) or spin co- 
herent state (SCS) |13j . rj is complex with unit modulus. 
From Eq. and the identity (-1)^1??) = | - ??), it is 
straightforward to show that the SU(2) ECS is the opti- 
mal input state which yields the maximal entanglement 
rate (3. 

The above analysis is restricted to finite equal- 
dimensional composite systems, but this limitation is 
convenient, not essential. Consider the harmonic oscilla- 
tor with infinite-dimensional Hilbert space. The operator 
a^a is the number operator, where (a) is the creation 
(annihilation) operator. We can have the Hamiltonians: 

H 2 = (-1)^® (-l) a+a and # 3 = (-l) at ° ® (-l) ot °. By 
noticing that the SCS can be realized in the Fock space 
as the binomial state 0], the optimal input states for 
Hamiltonians H^ and H% are directly obtained by ap- 
propriately replacing SCS with the binomial state in the 

SU(2) ECS CE9- 

Although in general one cannot expect an entangled 
state to be generated by Hamiltonian evolution from 
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a product state, in this spin system the optimal in- 
put state 1 77) ecs can be generated by the Hamiltonian 
—Hi from the product state \rj) (£> \rj). Mathematically 
|7?)ecs = exp(iHit)\ri} (g> |ry) with x = cos 2 t. This fact 
and the ancilla-independence of the self-inverse Hamilto- 
nian entanglement capability suggest that the entangle- 
ment capability depends more on the Hamiltonian and 
the associated unitary operator. Next, we demonstrate 
this character by introducing the entanglement rate for 
operator U(t), and give the result that the entanglement 
capability of Hamiltonian Hi is equal to the maximal op- 
erator entanglement rate for the associated unitary op- 
erator. 

We review operator entanglement introduced in 
Ref. 0. Notice that the linear operators over Tid (finite 
d) also form a c? 2 -dimensional Hilbert space denoted by 
Ti^p , and the corresponding scalar product between two 
operators X and Y is given by the Hilbcrt-Schmidt prod- 
uct (X,Y) := TrpT+Y), and ||A|| HS = (X,X). Then, the 
operator acting on Hd ± ®Hd 2 can be considered as a state 
on the composite Hilbert space H^ 2 S ® H^ 2 , and the op- 
erator entanglement can be defined as the entanglement 
of that state. 

Any unitary operator V acting on TLd t ® Ti.d 2 may be 
Schmidt-decomposed as 7} V = J2n s nA n ® B n , where 
s n > Vn and A n and B n are orthonormal operator 
bases for systems 1 and 2. From the Schmidt form the 
entanglement of a unitary operator V is determined to be 
£{V) — — J2i sf/{d\d2) log [sf /(did 2 )] , where the factor 
1/(^1^2) arises from normalization of the unitary oper- 
ator. We can think of the operator entanglement as a 
strength measure of the operator 0. 

If we consider the unitary operator U(t), the operator 
entanglement 8[U(t)] becomes a time-dependent func- 
tion. Analogous to the definition of the state entangle- 
ment rate |2|, it is natural to define the operator en- 
tanglement rate at a certain time of interaction and the 
maximal entanglement rate 

R(t) := d£[U(t)]/dt, R max := maxR(t), (19) 

respectively, where the maximization is over all time t. 
At certain times the entanglement rate becomes maxi- 
mal. 

The operator entanglement of an arbitrary unitary 
operator V is equal to the entanglement of the state 

£(V) = E[V(\$) A , A ®\$) B B>)}, (20) 

where 

n=Q 

\®)bb> = \ n )s ® \n)B> G H d2 ® H d2 (22) 

are maximally entangled states. Equation (|20|) shows a 
direct relation between operator entanglement and state 



entanglement. Note that here dimension di can differ 
from dimension g?2- From Eq. (|2U|) . we immediately have 
a relation between the operator entanglement rate and 
the state entanglement rate 

R(t) = r(t) = E[U(t)(\*) A >A ® \®)bb>)]- (23) 

From Eq. (|23|l we find that the maximal operator entan- 
glement rate gives a lower bound for the entanglement 
capability of a general Hamiltonian, which is the infinites- 
imal version of the one that the operator entanglement 
gives a lower bound of the entanglement capability of a 
unitary operator 0. 

Having defined the operator entanglement rate let us 
study the maximal operator entanglement rate i? max of 
the evolution operator associated with Hamiltonian Hi 
with dimension d = 2j + 1 an even number 15]. The 
associated unitary operator is in the Schmidt form 

Ui{t) = cosi l rf ®l rf -isint (-1)^ <g> (-1)^. (24) 

From the above equation, the operator entanglement and 
operator entanglement rate are given by 

E\Ui(t)]=- cos 2 t log cos 2 t - sin 2 t log sin 2 t, (25) 
R(t) =R[Ui{t)} = sin(2i) log(cot 2 t), (26) 

respectively. The above equation shows that the operator 
entanglement rate is a periodic function of time t; hence 
we can maximize over one period. It is straightforward 
to find that at time t — 0.2932, the entanglement rate 
reaches its maximum value, i.e., 

i?max - 13. (27) 

Comparing Eqs. (|16f) and (|27l) . we conclude that the en- 
tanglement capability of Hamiltonian Hi is equal to the 
maximum operator entanglement rate of the associated 
unitary operator exp(— iHit). Therefore, the entangle- 
ment capability of the self-inverse Hamiltonian is inher- 
ent in the evolution operator in the sense that it can be 
solely determined by the evolution operator, irrespective 
of states. 

More generally, there is interest in the entanglement 
capability of quantum o pera tors, not just Hamiltoni- 
ans 0, El [H QJ, [H IM M HI ■ Here, we analyze the 
entanglement capability of the unitary operator gener- 
ated by the self-inverse Hamiltonian H J2J • We quantify 
the entanglement capability of a unitary operator U by 
the maximum entanglement which a unitary operator can 
create given an initial product state 

E V = su P|7)i|5> £([/| 7 ) ® 1,5)), (28) 

where I7) G Ti A 'A and \5) € H. B B', namely, we include 
ancillas. For unitary operator U(t) we immediately ob- 
tain 

E u[t) =sup ljhls) E[U(t)\ 7 ) ® \5)} = 8up hm E[\9(t))], 

(29) 

|$(t)) =cosi \j)®\6) -isint \j)\S), (30) 
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with normalized states I7) = -JTaIt) and \5) = Xb\S). 

The state l|3Uf) belongs to a class of bipartite entangled 
states discussed in Ref. [23]. To quantify the entangle- 
ment of this state, we can use the standard entanglement 
measure, the entropy of entanglement. However, for con- 
venience, we use the concept of the concurrence C 24] to 
quantify the entanglement since the state can be viewed 
as a state of two pseudo-qubits |25j . Using the results of 
Ref. [25J, the concurrence of state {^(i)) is obtained as 

C[|*(t))] = I sin(2t)|^(l - |( 7 |7>l 2 )(l - 1ml 2 ) (31) 

From the above equation, it is easy to find that £r/(t) = 
|sin(2i)|, which is the maximal entanglement which the 
operator U{t) can generate from an arbitrary product 
state. As the states I7) e Ha'A and \5) e Hbb 1 , the 
entanglement capability Euu\ is ancilla-independent. 

The unitary operaror Ui(t) l|24l) is in fact a special case 
of U(t) J2J). Assuming that Ui(t) act on H d <S> H d with 
even d, we find that the optimal input state for generating 
maximal entanglement is \rf) ® \rf), and from Eq. (|24|) we 
find that 

£[U x (t)\ = C[U x {t)] = E Ul{t) = I sin(2t)|. (32) 



Therefore, the operator entanglement of U\ is equal to 
the entanglement capability of the unitary operator. 

In conclusion, we have determined the entanglement 
capability of self-inverse Hamiltonian evolution. The self- 
inverse Hamiltonians studied here are physically relevant 
in their own right, and reduce to the important Ising 
Hamiltonian that is central to quantum information sci- 
ence. These Hamiltonians go beyond the two-qubit cases; 
namely, they can act on d\ x c?2 composite system, where 
the dimension of the subsystem can be either finite or 
infinite. We introduced the concept of operator entangle- 
ment rate which is well-defined for composite finite sys- 
tems, and for certain self-inverse Hamiltonians the maxi- 
mal operator entanglement rate is equal to the entangle- 
ment capability. For a general Hamiltonian the maximal 
operator entanglement rate always gives a lower bound 
on the Hamiltonian entanglement capability. 
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